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Low mass lepton pair production in hadron collisions 
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On ', The hadroproduction of lepton pairs with mass Q and transverse momentum Qt can be described in pertur- 

^\ ■ bative QCD by the same partonic subprocesses as prompt photon production. We demonstrate that, like prompt 

On ' photon production, lepton pair production is dominated by quark-gluon scattering in the region Qt > Q/2. This 

leads to sensitivity to the gluon density in kinematical regimes that are accessible both at collider and fixed 

f-H ■ target experiments while eliminating the theoretical and experimental uncertainties present in prompt photon 

^ \ production. 
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1. INTRODUCTION 

The production of lepton pairs in hadron col- 
lisions hih2 -^ 7*X;7* — > ^f proceeds through 
an intermediate virtual photon and its subse- 
quent leptonic decay. Traditionally, interest in 
this Drell-Yan process has concentrated on lep- 
ton pairs with large mass Q which allows for the 
application of perturbative QCD and the extrac- 
tion of the antiquark density in the proton ||l| . 

Prompt photon production hih2 — *■ 'yX can 
be calculated in perturbative QCD if the trans- 
verse momentum Qt of the photon is sufficiently 
large. This process then provides essential in- 
formation on the gluon density in the proton 
at large x 0. Unfortunately, it suffers from 
considerable fragmentation, isolation, and intrin- 
sic transverse momentum uncertainties. Alterna- 
tively, the gluon density can be constrained from 
the production of jets with large transverse mo- 
mentum at hadron colliders H, which however 
suffers from ambiguous information coming from 
different experiments and colliders. 

In this paper we demonstrate that, like prompt 
photon production, lepton pair production is 
dominated by quark-gluon scattering in the re- 
gion Qt > Q/2. This leads to sensitivity to the 
gluon density in kinematical regimes that are ac- 
cessible both at collider and fixed target exper- 
iments while eliminating the theoretical and ex- 
perimental uncertainties. 



In Sec. H, we briefly discuss the relationship 
between virtual and real photon production in 
hadron collisions in next-to-leading order QCD. 
In Sec. |3| we present our numerical results, and 
Sec. y contains a summary. 



NEXT-TO-LEADING 
FORMALISM 



ORDER QCD 



* Supported by the U.S. Department of Energy, Division of 
High Energy Physics, under Contract W-31-109-ENG-38. 



In leading order (LO) QCD, two partonic sub- 
processes contribute to the production of virtual 
and real photons with non-zero transverse mo- 
mentum: qq — > 'y^*'g and qg — > 'y^*'q. The cross 
section for lepton pair production is related to 
the cross section for virtual photon production 
through the leptonic branching ratio of the vir- 
tual photon a/{3TTQ^). The virtual photon cross 
section reduces to the real photon cross section in 
the limit Q^ ^ 0. 

The next-to-leading order (NLO) QCD correc- 
tions arise from virtual one-loop diagrams inter- 
fering with the LO diagrams and from real emis- 
sion diagrams. At this order processes with inci- 
dent gluon pairs (gg), quark pairs (qq), and non- 
factorizable quark- antiquark (5^2) processes con- 
tribute also. Singular contributions are regulated 
in n = 4— 2e dimensions and removed through MS 
renormalization, factorization, or cancellation be- 
tween virtual and real contributions. An impor- 
tant difference between virtual and real photon 
production arises when a quark emits a coUinear 
photon. Whereas the coUinear emission of a real 
photon leads to a 1/e singularity that has to 



be factorized into a fragmentation function, the 
coUinear emission of a virtual photon gives a fi- 
nite logarithmic contribution since it is regulated 
naturally by the photon virtuality Q. In the limit 
Q^ — > the NLO virtual photon cross section re- 
duces to the real photon cross section if this loga- 
rithm is replaced by a 1/e pole. A more detailed 
discussion can be found in Q. 

The situation is completely analogous to hard 
photoproduction where the photon participates 
in the scattering in the initial state instead of the 
final state. For real photons, one encounters an 
initial-state singularity that is factorized into a 
photon structure function. For virtual photons, 
this singularity is replaced by a logarithmic de- 
pendence on the photon virtuality Q pi . 

3. NUMERICAL RESULTS 

In this section we present numerical results 
for the production of lepton pairs in pp colli- 
sions at the Tevatron with center-of mass energy 
^/S = 1.8 and 2.0 TeV and in pH"^ coUisions at 
fixed target experiments with \/S = 38.8 GeV. 
We analyze the invariant cross section Ecfia/dp^ 
averaged over the rapidity interval -1.0 < y < 1.0 
at the Tevatron and averaged over the scaled lon- 
gitudinal momentum interval 0.1 < xp < 0.3 at 
fixed target experiments. We integrate the cross 
section over various intervals of Q and plot it as 
a function of the transverse momentum Qt- Our 
predictions are based on a NLO QCD calcula- 
tion and are evaluated in the MS renormaliza- 
tion scheme. The renormalization and factoriza- 
tion scales are set to fi = iif = \J Q"^ + Q\. If 
not stated otherwise, we use the CTEQ4M par- 
ton distributions pX| ] and the corresponding value 
of A in the two-loop expression of ois with four 
fiavors (five if /x > m(,). The Drell-Yan factor 
Q;/(37r(3^) for the decay of the virtual photon into 
a lepton pair is included in all numerical results. 

In Fig. |l| we display the NLO QCD cross sec- 
tion for lepton pair production at the Tevatron 
at \fS = 1.8 TeV as a function of Qt for four 
regions of Q. The regions of Q have been cho- 
sen carefully to avoid resonances, i.e. between 
the p and the J li\j resonances, between the J/V' 
and the T resonances, above the T's, and a high 



pp^y'XalVs= 1.8 TeV 



2 GeV < Q < 3 GeV 



5 GeV < Q < 6 GeV 



11 GeV<Q< 12 GeV 



30 GeV < Q < 35 GeV 



-1 <y< 1 



10 20 30 40 50 60 70 80 90 100 
Qt (OcV) 

Figure 1. Invariant cross section E(Pa/dp^ as a 
function of Qt for pp -^ j*X at ^/S = 1.8 TeV 
in non-resonance regions of Q. The cross section 
falls with the mass of the lepton pair Q and, more 
steeply, with its transverse momentum Qt- 



mass region. The cross section falls both with the 
mass of the lepton pair Q and, more steeply, with 
its transverse momentum Qt- Unfortunately, no 
data are available yet from the CDF and DO ex- 
periments. However, data exist for prompt pho- 
ton production out to Qt — 100 GeV, where the 
cross section is about 10~^ pb/GeV^. It should 
therefore be possible to analyze Run I data for 
lepton pair production up to at least Qt — 30 
GeV where one can probe the parton densities in 
the proton up to xt = 2Qt/VS ~ 0.03. The 
UAl collaboration measured the transverse mo- 
mentum distribution of lepton pairs at VS — 630 



GeV up to XT = 0.13 |11|, and their data agree 
well with our theoretical results Q . 

The fractional contributions from the qg and 
qq subprocesses up through NLO are shown in 
Fig. 0. It is evident from Fig. that the qg sub- 



reactions begins to raise the fraction of the cross 
section attributed to the qq subprocesses. 

Data obtained by the Fermilab E772 collabo- 
ration [|2| from an 800 GeV proton beam inci- 
dent on a deuterium target are shown in Fig. ^ 
along with theoretical calculations. For our anal- 
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Figure 2. Contributions from the partonic sub- 
processes qg and qq to the invariant cross section 
E(P(j/dp^ as a function of Qt for pp -^ 'y*X at 
^/S = 1.8 TeV. The qg channel clearly dominates 
in the region Qt > Q/2. 



process is the most important subprocess as long 
as Qt > Q/2. The dominance of the qg subpro- 
cess diminishes somewhat with Q, dropping from 
over 80 % for the lowest values of Q to about 70 
% at its maximum for Q — 30 GeV. In addition, 
for very large Qt, the significant luminosity asso- 
ciated with the valence dominated q density in pp 




Figure 3. Invariant cross section Ed?a/dp^ as a 
function of Qt for pN ^-f*X at pub = 800 GeV 
in the region between the J/ip and T resonances. 
The NLO perturbative cross section (solid) is 
shown along with the all-orders resummed expec- 
tation (dashed), the asymptotic result (dotted), 
and a matched expression (dot-dashed). The data 
are from the Fermilab E772 collaboration. 



ysis we have chosen a lepton pair mass region 
between the J/tp and T resonances. The solid 



line shows the purely perturbative NLO expecta- 
tion. The transition to low Qt can be described 
by the soft-gluon resummation formalisin and is 
shown in the dashed curve in Fig. ^ |QJ8) . The re- 
summed result can be expanded in a power series 
in as asymptotically around Qt = 0. Its NLO 
component (dotted curve) can then be matched 
to the perturbative result (dot-dashed curve) ||]. 
From Fig. || it becomes clear that resummation 
is not needed and fixed order perturbation the- 
ory can be trusted when Qt > Q/2. Unfortu- 
nately, the data from E772 do not extend into 
this region. However, the cross section should be 
measurable in forthcoming experiments down to 
10-3 pb/GeV^, i.e. out to at least Qt ^ 6 GeV 
or XT — 0.31 where the gluon density is poorly 
constrained now. 

In Fig. H we demonstrate that in fixed target 
experiments also the lepton pair cross section is 
dominated by quark-gluon scattering at the level 
of 80 % once Qt — Q- The results in Fig. | also 
prove that subprocesses other than those initiated 
by the qq and qg initial channels are of negligible 
import. 

We will now turn to a previously unpublished 
study of the sensitivity of collider and fixed target 
experiments to the gluon density in the proton. 
The full uncertainty in the gluon density is not 
known. Here we estimate this uncertainty from 
the variation of different recent paramctrizations. 
We choose the latest global fit by the CTEQ col- 
laboration (5M) as our point of reference Pi and 
compare it to their preceding analysis (4M px|]) 
and to a fit with a higher gluon density (5HJ) in- 
tended to describe the CDF (and DO) jet data 
at large transverse momentum. We also com- 
pare to global fits by MRST [^, who provide 
three different sets with a central, higher, and 
lower gluon density, and to GRV98 |l|]0 For this 
study we update the Tevatron ccnter-of-mass en- 
ergy to Run II conditions {\/S = 2.0 TeV) which 
increases the invariant cross section for the pro- 
duction of lepton pairs with mass 5 GeV < Q < 



■^In this set a purely perturbative generation of heavy 
flavors (charm and bottom) is assumed. Since we are 
working in a massless approach, we resort to the GRV92 
parametrization for the charm contribution [h4| and as- 
sume the bottom contribution to be negligible. 



6 GeV by 5 % at low Qt :^ 1 GeV and 20 % at 
high Qt =i 100 GeV. 

In Fig. ^ we plot the cross section for lepton 
pairs between the J/tp and T resonances at Run 
II of the Tevatron which should be measurable up 
to at least Qt =i 30 GeV (xt ^ 0.03). For the 
CTEQ paramctrizations we find that the cross 
section increases from 4M to 5M by 2.5 % {Qt = 
30 GeV) to 5 % {Qt = 10 GeV) and from 5M to 
5HJ by 1 % in the whole Qr-range. The largest 
differences to CTEQ5M are obtained with GRV98 
at low Qt (minus 10 %) and with MRST(gt) at 
large Qt (minus 7%). 

A similar analysis for conditions as in Fermi- 
lab's E772 experiment is shown in Fig. H. In 
fixed target experiments one probes substantially 
larger regions of Xt than in collider experiments. 
Therefore one expects a much larger sensitivity to 
the gluon distribution in the proton. Indeed we 
find that CTEQ5HJ increases the cross section 
by 7 % (26 %) w.r.t. CTEQ5M at Qt = 3 GeV 
{Qt = 6 GeV) and even by 134 % at Qt = 10 
GeV. For MRST(gi) the CTEQ5M cross section 
drops by 17 %, 40 %, and 59 % at these three 
values of Qt ■ 



4. SUMMARY 

In summary, we have demonstrated that the 
production of Drell-Yan pairs with low mass 
and large transverse momentum is dominated 
by gluon initiated subprocesses. In contrast to 
prompt photon production, uncertainties coming 
from fragmentation, isolation, and intrinsic trans- 
verse momentum are absent. The hadroproduc- 
tion of low mass lepton pairs is therefore an ad- 
vantageous source of information on the gluon 
density in the proton at large x in collider experi- 
ments and even more in fixed target experiments. 
Massive lepton pair production data could pro- 
vide new insights into the parametrization and 
size of the gluon density. 
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Figure 4. Contributions from the NLO QCD par- 
tonic subprocesses to the invariant cross section 
Ed^a/dp^ as a function of Qt for pN —> j*X 
at piab — 800 GeV. qg (solid) dominates over qq 
(dashed) at the level of 80 % once Qt — Q, and 
the pure NLO QCD processes gg (dotted), qq2 
non-factorizable parts (dot-dashed) , and qq (wide 
dots) are negligible. 
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Figure 5. Invariant cross section Eifia/dp^ as 
a function of Qt for pp -^ ^*X at ^fS = 2.0 
TeV in the region between the J /tp and T res- 
onances. The largest differences from CTEQ5M 
are obtained with GRV98 at low Qt (minus 10 %) 
and with MRST(gt) at large Qt (minus 7 %). 



pN^YXatP|^j = 800GeV 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


'%■ 


>■'%, CTEQ5M 


^ ' \^-. 


'\''-\ 


^'■\. 


^.■-%, CTEQ4M 




\ ■ M'' rf^'I'T 'f^CJ TT 


_ ■ V V, LlhQ5HJ_ 


\ '"x\'-. 


"~'. "^^■■■, 


\^\\, MRST - 


\ ~\\''-. 


\, ''X^'-.. 


\'\\\ MRST(gt)" 




\ ~- x''-. 


\ '\\\ MRST(g4.) 


\''% 


\ ~-. \N- GRV98 


^. ■-■>■■, 


\ -\ ■■■ 


\ ■■%■■■■■■ 


^. \\\ 


\ ' '.\ ''■■■■ 




'^ '~\\-'-'-. 


\ \ V '■., 


~ \ \ \,''-- ~ 


■^ '- \- '■■. 


\, "■, \-''-., 


5GeV<0<6GeV ^^ '^\\ " 


'\ ~\ \' '"\ 


''■ '■' \ '"■■ 


0.1<Xp<0.3 \ ■\ \ 


\ ^^ ^ 


^, 



3 3.5 4 4.5 5 5.5 6 

Qt (GcV) 

Figure 6. Invariant cross section E(Pa/dp^ as a 
function of Qt for pN ^-f*X at piab = 800 GeV. 
The cross section is highly sensitive to the gluon 
distribution in the proton in regions of xt where 
it is poorly constrained. 



